Boost-invariant mean field approximation and the nuclear Landau-Zener effect 
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We investigate the relation between time-dependent Hartree-Fock (TDHF) states and the adi- 
abatic eigenstates by constructing a boost-invariant single-particle Hamiltonian. The method is 
numerically realized within a full three-dimensional TDHF which includes all the terms of the 
Skyrme energy functional and without any symmetry restrictions. The study of free translational 
motion of a nucleus demonstrates the validity of the concept of boost-invariant and adiabatic TDHF 
states. The interpretation is further corroborated by the test case of fusion of 16 0+ 16 0. As a first 
application, we present a study of the nuclear Landau-Zener effect on a collision of 4 He+ ls O. 
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INTRODUCTION 



The time-dependent Hartree-Fock (TDHF) approx- 
imation, originally proposed by Dirac [H, has found 
widespread applications in nuclear dynamics since more 
than thirty years 0]. It provides the microscopic foun- 
dation for describing various dynamical scenarios in the 
regime of large amplitude collective motion. Soon after 
its introduction to nuclear physics, the TDHF approxi- 
mation was extensively applied to studies of fusion exci- 
tation functions, fission, deep-inelastic scattering, collec- 
tive excitations and nuclear molecular resonances. These 
studies shed light on several mechanisms of heavy-ion col- 
lisions and resonance dynamics. Reviews on these earlier 
TDHF applications can be found, e. g., in Ref. 0,|3|- 

At that time, however, limited computer capacity re- 
stricted most calculations to axial symmetry and omis- 
sion of spin-orbit coupling. These limitations turned 
out to be a hindrance for the development. For exam- 
ple, earlier TDHF calculations underestimated the en- 
ergy dissipation from the collective kinetic energy into 
internal excitations so that the energy window of fu- 
sion reactions was too small in comparison with experi- 
ments. Later work demonstrated that the spin-orbit cou- 
pling 0, Q and fully three-dimensional geometry @, H, Q 
enhance dissipation. None of these calculations, how- 
ever, was able to include all constituents simultane- 
ously. With the steady upgrade of computational power, 
three-dimensional TDHF calculations employing the full 
Skyrme force became possible and renewed the interest 
in nuclear TDHF as seen from recent publications on res- 
onance dynamics fToL [Til [l2l . [Til ITTj and heavy-ion colli- 
sions [l5|, [l6|, [Tt], Il8j . One expects that the new genera- 
tion TDHF calculations may yield more realistic features 
for heavy-ion collisions at low and medium energies and 
for resonance dynamics. This revives old questions that 
have been left unanswered for a while, one of which is the 
subject of this paper: the analysis of heavy-ion collisions 
as computed by TDHF in terms of adiabatic states and 
level-crossing dynamics. 

Nucleus-nucleus collisions present different behavior 



depending on the delicate balance of reaction time and 
rearrangement time of the mean field. For example, the 
experimentally observed resonance-like peaks in the in- 
clastic cross sections are explained in terms of the well- 
known Landau-Zener excitation due to a breakdown of 
the adiabaticity condition near an avoided level cross- 
ing. The Landau-Zener mechanism was first introduced 
into nuclear physics in Ref. [l~9j . The Landau-Zener ef- 
fect and its a pplic ations to heavy-ion collisions have been 
discussed pfj l2ll [22j in terms of the asymmetric two- 
center shell model (TCSM) [23[, employing the assump- 
tion that the nucleons can be described with adiabatic 
"molecular" states, i. e. the instantaneous eigenstates of 
the deformed mean field, during the heavy ion reaction. 
A review on Landau-Zener dynamics and experimental 
data in nuclear molecules may be found in Refs. [24ll25l|. 

TDHF calculations, however, have never been analyzed 
in terms of the Landau-Zener effect because their single- 
particle states and energies have no simple physical inter- 
pretation and a construction of corresponding adiabatic 
states was not readily available. There were early at- 
tempts to define a related adiabatic basis by means of 
density-constrained Hartree-Fock [26|, [2?j] ■ These turned 
out to be very promising, lacking, however, at that time 
the exact treatment of the flow contributions. Improved 
computing power now allows revisiting the case with- 
out technical restrictions. We thus will present here a 
self-consistent scheme to define and compute two useful 
analyzing instruments for a given TDHF state: instan- 
taneous single-particle energies and instantaneous adia- 
batic states. The straightforward expectation values of 
the single-particle Hamiltonian turn out to be blurred by 
trivial flow contributions. To take into account the effect 
of motion on the single-particle wave functions, we define 
a single-particle Hamiltonian which is invariant under 
Galilei transformation, in particular under a boost. We 
call that the boost-invariant Hamiltonian. This dramat- 
ically reduces the energy variances of the actual TDHF 
states, providing better-defined single-particle energies. 
Moreover, the (instantaneous) eigen-states and energies 
of the boost-invariant Hamiltonian provide a well de- 
fined adiabatic basis. For example, the single-particle 
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states in free translational motion are exact eigenstates 
of this boost-invariant Hamiltonian and its expectation 
values remain the static single-particle energies. While 
this property is not exact anymore in a situation of two 
colliding nuclei, we shall demonstrate that this method 
nevertheless allows to define meaningful single-particle 
energies and variances thereof. Moreover, one can estab- 
lish an approximate relation of the time-dependent solu- 
tions to the adiabatic deformation-dependent spectrum. 
As a first application we study the Landau-Zener effect 
in heavy-ion collisions. 

The paper is organized as follows: Section |TT] briefly 
recalls the Skyrme energy functional and TDHF. In sec- 
tion Hn] we construct a boost-invariant Hamiltonian and 
test its validity for free translational motion of a nucleus. 
Section IIVI applies the newly developed scheme to the 
analysis of level crossings in heavy-ion collisions. Section 
IVlis devoted to the summary and conclusion. 

II. COMMENTS ON TDHF WITH SKYRME 
FORCES 

Most nuclear TDHF calculations by far are based on 
the Skyrme energy functional; for a recent extensive re- 
view see [28| . It is also used in the present applica- 
tions. The starting point is the Skyrme energy-density 
functional Ssk = £{p, T i j\ J), which is expressed in 
terms of a few local densities and currents: density p, 
kinetic density r, spin density a, current j, and spin- 
orbit density J. It includes free kinetic energy, Skyrme 
interaction, Coulomb energy, and the center-of-mass cor- 
rection. The pairing energy is ignored here as we will 
deal with collisions of closed-shell nuclei. There are var- 
ious parametrizations of the Skyrme force (28j . Since 
the present study is concerned with fundamental effects 
which should not depend on the detailed force used, we 
thus chose just one parametrization out of many, namely 
the force SLy6 [29[ which is widely used and provides a 
reliable description of nuclear structure and dynamics. 

Using the principle of least action and varying with re- 
spect to the single-particle state <p* , we obtain the TDHF 
equations (in the following units with K = 1 are used) 

id t ip a = h(p,T,a,j,J)tp a , (1) 

where h is the time-dependent mean-field Hamiltonian 
depending on the occupied single-particle wave functions 
through densities and currents. Given the initial condi- 
tions, {ip a (r, t = 0)}, the TDHF equations (JTJ determine 
the wave functions for all later times. In the stationary 
limit, we obtain the static mean- field equation 

hip a = £ a lfi a , (2) 

where the single-particle energies e a appear naturally as 
eigenvalues of the mean-field Hamiltonian h. The ques- 
tion is how to generalize the definition of the single- 
particle energy to TDHF. The naive definition is to use 



the expectation value of the instantaneous mean field 
h(t). This, however, raises difficulties as we will see. Pos- 
sible improvements will be developed in the sequel. 

A few words on the numerical solution are in order. 
The set of non-linear TDHF equations is solved on a 
three-dimensional Cartesian coordinate-space grid em- 
ploying a Fourier representation for the derivatives. All 
contributions of the full Skyrme force were included and 
no symmetry restrictions imposed. The coordinate-space 
grid consists of 24 x 24 x 24 points with a grid spacing of 
1 fm. For the dynamical time stepping, we use a Taylor 
series expansion of the unitary mean-field propagator up 
to sixth order [f| and a time step of 0.2 fm/c. These nu- 
merical parameters provide good conservation of particle 
number and total energy during the dynamic evolution. 
The static HF equations were solved with the damped 
gradient iteration method [3(], HH . 



III. SINGLE-PARTICLE ENERGIES IN A 
MOVING FRAME 

A. Adiabatic expansion as a propaedeutic example 

Adiabatic single-particle states arc eigenstates of a 
single-particle Hamiltonian for a given set of deforma- 
tion parameters. For a first introduction, we will discuss 
that concept in this section on the grounds of a given, 
properly parametrized, sin gle- particle Hamiltonian, such 
as, e. g., that of the TCSM [23|]. Such a Hamiltonian usu- 
ally depends on a few collective deformation parameters 
which characterize the wanted reaction path. For sim- 
plicity let us just deal with internuclcar distance R as the 
sole such parameter and skip trivial complications such 
as spin and isospin. The proper choice of an independent 
single-particle potential V(f; R) and the subsequent so- 
lution of the eigenvalue problem yield a set of adiabatic 
single-particle states 4>k {f\ R) and corresponding eigenen- 
ergies €k(R). Using the adiabatic single-particle states 
the time-dependent (but still independent-particle) solu- 
tion is expanded as 

^(r, t) = J2 c jk (t)Mr; R{t))e- { ti dt ' ^ M . (3) 
k 

Such an expansion underlies, e. g., the cranking 
model (321 . I33T ]. This expansion has the problem that 
it relies on a stationary basis in which the current van- 
ishes for all states. Any flow has to be described through 
the complex expansion coefficients Cjk- This limits the 
ansatz ([3]) to extremely slow motion. 

An instructive example is uniform center-of-mass 
translation of an unexcited nucleus with velocity v = 
p/m. It is described by coherently moving single particle 

states ijjk(r, t) = 4>k(r — vt) cxp (i(k ■ r— efet) J . Clearly, 

the trivial plane-wave factor which produces the cor- 
rect flow is missing from the basis states in It has 
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to be reconstructed laboriously by the expansion coeffi- 
cients. A much more efficient description could be ob- 
tained by properly extending the scheme to a dynamic 
basis (j>k (r; R, R) which already accounts for collective 
flow. This step was found to be crucial for the deriva- 
tion of microscopic theories for collective motion in the 
framework of adiabatic TDHF [34| . We will now use the 
extension for the definition of adiabatic reference states 
in TDHF. The example of center-of-mass motion will be 
used as guidance. 



h. The variance grows quadratically with the boost mo- 
mentum P, i. e., proportionally to the center-of-mass en- 
ergy. The expectation value i a also becomes misleading. 
The kinetic contribution makes the binding properties 
invisible. This problem was already noticed by Thou- 
less and Valatin [35| while they were studying Galilean 
invariance of the TDHF equation. 

C. Construction of a boost-invariant mean field 



B. Flow-induced variances 

We will now discuss the case of self-consistent mean 
fields. To simplify the formal considerations, we restrict 
the discussion to one spatial dimension and think in 
terms of the simplest energy functional £ = £(p) depend- 
ing only on the local density and producing a purely local 
mean field, i. e., 

v 2 5£ 
h[p] = ^ + U(x,t) , U( X ,t) = WJ) , (4) 

where U is obtained from £ by functional derivation. The 
ground-state wave functions {<^o,oj fulfill the equations 

(ho - e a )<Po,a(x) = (5) 

with e a being the static single-particle energies. 

As a test case, consider center-of-mass motion of the 
HF ground state. The motion is initialized by a boost 
with total momentum P = MV where V is the velocity 
of the center of mass and M = Nm the total mass. The 
same boost is applied to all single-particle wave functions 

^o,a — > <p a (x,t) = e lP& l N y , a (x - Vt)e-^\ (6a) 

p2 

ia = £a+ 2M ■ (6b) 

The local density is then propagated with velocity V as 
p(x,t) = po(x — Vt) and this, in turn, carries through 
to the mean field motion as U(x,t) = Uo(x — Vt). 
The boosted wave functions together with the similarly 
moving mean field are the solution of the TDHF equa- 
tions (JTJ). The action of the mean field on the boosted 
wave function can be expressed in terms of the static 
solution as 

hp a = e' pi l N e-^ t (e a + ^- i p)^ a . (7) 

The expectation value is simply ((p a \h\ip a ) = e a . Both 
together allow evaluation of the variance of the single- 
particle Hamiltonian explicitly as 

P 2 

(tp a \Ah 2 \tp a ) = —(ip 0a \p 2 \ip 0a ). (8) 

It is obvious that the moving wave function ip a is not an 
cigenstate of the instantaneous mean-field Hamiltonian 



The above example of center-of-mass motion is instruc- 
tive. The variance of the mean-field Hamiltonian grows 
although we know that the system remains intrinsically 
unaltered. All that happens is a trivial kincmatical effect. 
Thus there should be ways to undo it equally trivially. 
In the center-of-mass case, we could simply transform 
the single-particle momentum into the intrinsic frame as 
p — ► p — P/N and use that in the kinetic energy oper- 
ator. That indeed provides a reasonable boost-invariant 
Hamiltonian for that particular case. A generalization 
can be obtained with the concept of the local momen- 
tum distribution p(x) as given by the local current j(x). 
This suggests the definition of a locally boost-invariant 
momentum 

J (x) 

P ► Pinv =P~ P~(x) , p(x) = — , (9) 

p(x) 

which can be extended to a boost-invariant kinetic- 
energy density r inv = J2 a |Pinv<p Q | 2 = r - j 2 /p. It is 
interesting to note that this is practically the Galilean- 
invariant combination rp — j 2 of kinetic contributions 
in the interaction part of the Skyrme energy func- 
tional [28|, [36| . This gives confidence in the above gener- 
alization. 

The idea thus is to define an "intrinsic" energy- 
functional by replacing the kinetic energy oc J d 3 rr by 
the boost-invariant kinetic energy 

£ kin , inv = ^J dx ( T - J —) ■ ( 10 ) 

The potential energy was already boost-invariant and 
thus the total functional becomes so. This functional 
is to be used for the purpose of analysis only, and it 
plays no role for the computation of the time evolution 
as such. Variation leads to the corresponding boost- 
invariant mean-field Hamiltonian 

where {...,...} is the anti-commutator and U(x,t) the 
usual time-dependent mean-field potential. The first two 
terms are exactly the same as in the usual TDHF Hamil- 
tonian and the last two stem from the boost-invariant 
kinetic energy. The corrected Hamiltonian then defines 
a boost-invariant single-particle energy 

4 inv) = (<p a \h inv \tp a ). (12) 
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Next we will show that the boost-invariant Hamiltonian 
has the boosted TDHF wave functions as eigenstates in 
the case of free translation. 

D. Test case: free translational motion 

The construction of the boost-invariant Hamiltonian 
(fTTj) was guided by the example of free center-of-mass 
translation, so that the natural test case is the global 
center-of-mass boost (pa|) . In that case, the local mo- 
mentum distribution is given by 

_ j(x) P 

PO) = —TT = T7 = const. , (13) 
pyx) Jy 

while the boost-invariant Hamiltonian reduces to 

p 2 P P 2 

h ™ = £- + U(x,t)-—p+— , (14) 
2m M 2M 

and finally the expectation value of the boost-invariant 
Hamiltonian and its variance become 

(ip a \h inv \ip a ) = e a , (15a) 
(<p a \Ah 2 nv \<p a ) = 0, (15b) 

where e a are the static single-particle energies. 

The zero variance means that the boosted TDHF wave 
function ip a according to Eq. (|6ap is an eigenstate of 
the boost-invariant Hamiltonian h- lnv and its eigenvalue 
remains static solution e a as defined in Eq. ([5]). 

These results for free translation suggest the boost- 
invariant Hamiltonian h- lm as an appropriate instrument 
for analyzing the single-particle states of TDHF. The 
single-particle energies computed as expectation values 
of the boost-invariant Hamiltonian (jlip can be consid- 
ered as "intrinsic" single-particle energies representing 
the actual binding independent from trivial kincmatical 
contributions. It is to be noted that these equations are 
also applicable for a non-local mean-field Hamiltonian 
like, e.g., for the Skyrme force. 

The practical computation of the boost-invariant 
Hamiltonian is a bit demanding due to the density in 
the denominator. Nevertheless, for free translation of 
the nucleus 16 we achieve variances of about 0.02-0.05 
MeV and expectation values stay within 10 -4 ~10 -5 MeV 
of the static ones with the full three-dimensional TDHF 
which includes all the terms of the Skyrme energy func- 
tional and without any symmetry restrictions. 
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FIG. 1: Time evolution of boost-invariant and adiabatic ob- 
servables for a central 16 0+ ie O collision at 25 MeV. The 
dashed lines with error bars show the boost invariant single- 
particle energies ei mv ' and their variances Ael mv ' for the low- 
est (lsi/2) proton states, and solid lines the adiabatic energies 

(adia) r 

e\ for comparison. 

then become something different. They correspond to 
some extent to the adiabatic states and we will therefore 
use the expressions "adiabatic states" and "adiabatic en- 
ergies" in the following. A given TDHF state ip a imbued 
with some intrinsic excitation is distributed over the adi- 
abatic basis {4>i}- This can be quantified through the 
adiabatic occupation probability 

^ (OCC) = E KVa|fc>| 2 (17) 
aEocc 

where the sum runs over the occupied TDHF states tp a . 
That quantity is the probability to find the state <f>i oc- 
cupied when expanding the actual TDHF Slater state 
into the adiabatic basis. Complementarily we have the 
hole probability p( unocc '> = 1 _ p( occ \ ^he occupation 
probability quantifies in its way the amount of intrinsic 
excitation carried in the TDHF states. The value of one 
means no excitation at all for states of the adiabatic basis 
below the Fermi level, and lowering below one is closely 
related to excitation. 



IV. ANALYSIS OF HEAVY-ION COLLISIONS 

A. Level schemes 



E. Adiabatic states 

True intrinsic excitations in more general dynami- 
cal situations add some energy variance to the TDHF 
states. The eigenstates of the instantaneous boost- 
invariant Hamiltonian, 

hn^i = <h 4>i , (16) 



For the investigation of heavy-ion collisions, the wave 
functions of the two fragments are placed symmetrically 
on the grid 5 fm off the center of box and then boosted to 
the desired relative center-of-mass energies. This is the 
initial state for TDHF dynamical propagation. Besides 
the calculation of the usual TDHF states {(p a } of the 
colliding system, diagonalization of the boost-invariant 
Hamiltonian gives the adiabatic single-particle energies 
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FIG. 2: Time evolution of higher lying proton states for a central collision of O+ O with a center of mass energy of 25 MeV. 
The upper panels show the boost invariant single-particle energies ei mv ' (dashed) and the adiabatic energies e ' adla ) (f u ll lines). 
The lower panels show the corresponding adiabatic occupation probabilities. The arrows pointing from the upper panels to 
the lower ones indicate the level crossings where occupation changes rapidly. The labels of the single-particle states at t=0 are 
indicated near the left axis. 



and wave functions {4>i}- It should be noted that in a col- 
lision situation the boost-invariant Hamiltonian cannot 
be expected to eliminate the kinetic effects completely, 
because it depends on the local conditions. It is hoped, 
however, that for the initial stages of the interaction the 
behavior of the levels can still be extracted. For the sym- 
metric system, parity projection has been done for the 
boost-invariant wavefunctions. 

To illustrate the performance of boost-invariant and 
adiabatic states, Fig. [1] shows a typical result for the 
splitting of an initially highly degenerate state during 
a heavy-ion reaction. The test case is a central 16 0+ 16 
collision with center-of-mass energy of 25 MeV. Shown 
are the two lowest levels (the proton ls 1 / 2 levels in the 
right and left collision partner). The first few fm/c show 
the initial phase with nearly free cm. motion of the two 
nuclei towards each other. Boost-invariant and adiabatic 
states remain identical and the energy variance is zero 
within the limit of numerical precision. The nuclei start 
to interact around 30 fm/c. At this time the initially 
degenerate lsj/2 states split into two levels, as expected. 
The way the levels split into the various sub-states is 
very similar to what is expected from a two-center ap- 
proach like in [23| . At the same time the energy variance 
of the boost-invariant states grows because the collision 
mixes forward with backward flow such that we have an 
increasing spread of flow around its decreasing average. 
It is to be noted that the variances of the boost-invariant 
Hamiltonian are much smaller than those of the usual 
TDHF Hamiltonian. The latter comes up to several tens 



of MeV clearly showing that the single-particle energies 
from the TDHF states have no meaning at all. Further- 
more, the boost-invariant energies become slightly larger 
than the adiabatic ones, which expresses the amount of 
true intrinsic excitation piling up in the boost-invariant 
states. It is very satisfying to see that the energy differ- 
ence is proportional to the variance of the boost-invariant 
states. The signals confirm each other as measure for in- 
trinsic excitation and they give credibility to both forms 
of energy expectation values. Altogether, figure Q] clearly 
demonstrates the usefulness of the concept and the rela- 
tion between the boost-invariant states containing local 
flow and the quasi-stationary adiabatic states. 

Figure [2] shows the time evolution of the other occu- 
pied proton states in our example of a 16 0+ 16 collision. 
As for the lsi/ 2 levels above, we again see the splitting 
of the asymptotically degenerate states \pz/2 an d lPi/2 
with increasing interaction. The difference between adi- 
abatic and boost-invariant energies indicates the degree 
of internal excitation. It can be very different for the 
different states. It is, e. g., somewhat surprising that 
the lower state (left panels) acquires excitation rather 
early while some higher states wait for much longer. The 
lower panels of figure[2]show the occupation probabilities 
(jTTJ). They start at unity as they should for a yet un- 
excited state where adiabatic and boost-invariant states 
are still identical. Their subsequent decrease reflects the 
degree of intrinsic excitation. The two indicators for in- 
trinsic excitation, occupation probability and difference 
e (mv) _ £ (adia)^ a g ree n j ce iy f or & \\ states shown. The adi- 
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abatic energies in the upper panels of figure [2] show the 
interesting phenomenon of level crossings, quite similar 
to that observed in former studies using deformed shell 
models [20|, HH HH . The arrows connect these points 
with the lower panels where the occupation probabilities 
seem to jump. But that is merely a labeling effect when 
sorting the states always according to adiabatic energy. 
Following the states diabatically through the crossings 
would produce smooth evolution of energies and occupa- 
tion probabilities. Such a diabatic tracking, however, is 
only possible if we ignore pairing, as we do here. Inclu- 
sion of pairing would smoothen the crossings and enforce 
adiabatic tracking with subsequent "smooth jumps" in 
the observables. 



B. Landau-Zener effect 

As mentioned in the introduction, reactions of com- 
plex many-body systems like nuclei or molecules often 
produce the Landau-Zener effect. It happens at level 
crossings with only small coupling between the two en- 
ergetically close levels. There is a competition between 
the speed with which the levels evolve and the time nec- 
essary for the rearrangement of the occupations. Very 
slow motion leaves sufficient time such that always the 
energetically lower level is fully occupied. That is the 
adiabatic limit. Increased velocity causes transitions (= 
diabaticity) where the occupation partially crosses over 
into the then higher level, thus turning collective en- 
ergy to internal excitation. That is the much celebrated 
Landau-Zener effect. In this section, we will employ the 
boost-invariant and adiabatic states as analyzing tools 
to a study of the Landau-Zener effect in a self-consistent 
mean-field description of heavy-ion collisions. For this 
purpose an asymmetric reaction is more appropriate and 
we select a head-on collision of 4 He+ 16 0. 

Fig. 03a) shows the root-mean-square (RMS) radius of 
the colliding system as a function of time. It is clear 
that both the collisions with 75 and 125 MeV are deep- 
inelastic scattering. The adiabatic occupation probabil- 
ity as defined in Eg. (TIT)) for the last occupied neutron 
state labelled at initial stage as lpi/2 is presented in 
Fig-EIb). The assignment "occ." and "unocc." in Fig. [3] 
refers to the situation in the initial stage, while in the col- 
liding region the adiabatic occupation probability gives 
information on the actual occupation. 

The occupation probability is, of course, nearly unity 
during the initial stage of nearly free translation and 
starts decreasing as the two colliding nuclei approach 
each other. It then increases again and finally returns 
to nearly unity, not quite attaining it owing to the small 
transfer and evaporation probabilities after the separa- 
tion. The corresponding adiabatic single-particle ener- 
gies of the highest-lying occupied and the lowest-lying 
unoccupied neutron states, labelled at the initial stage 
as lpi/2 and ld 5 / 2 respectively, are shown in Fig. EHc). 
In the initial stage of dynamic propagation, the adiabatic 



single-particle energies are almost the same as those of 
the static ground states (realized exactly at t = 0). 

This behavior is easy to understand since the boost- 
invariant and adiabatic states only reflect the excitation 
and interaction of colliding nuclei. The two adiabatic 
states display the feature of avoided crossing around the 
smallest distance of two colliding nuclei. The same fea- 
ture also appears in the adiabatic occupation probabil- 
ity. The right panel with larger bombarding energy of 
125 MeV shows that the mixing of occupied and unoc- 
cupied components in the colliding stage becomes much 
stronger such that the occupied and unoccupied states 
are exchanged completely. Since the two interacting 
single-particle states belong to the nucleus 16 O, we find 
that two neutrons are excited from the uppermost oc- 
cupied state to the lowest unoccupied state. This ex- 
citation is activated gradually with increasing incident 
energy. This is a clear signal of a nuclear Landau-Zener 
transition in the TDHF description of a deep inelastic 
collision. 



V. SUMMARY 

In this work, we have constructed a boost-invariant 
single-particle Hamiltonian to eliminate the dynamically 
induced variances coming from the local velocity field in 
TDHF. For the case of free translational motion of a nu- 
cleus, the boost-invariant Hamiltonian produces eigen- 
states which have zero dynamical variances and repro- 
duce the stationary single-particle energies. In the case 
of a reaction, true intrinsic excitations take place and 
the TDHF states do not remain eigenstates of the boost- 
invariant Hamiltonian anymore. Their variances then be- 
come a measure of intrinsic excitation and the expecta- 
tion values still remain useful measures of single particle 
energies. Moreover, the eigenstates of the boost-invariant 
Hamiltonian can be considered as the (instantaneous) 
adiabatic states which contain no flow. The relation be- 
tween boost-invariant and adiabatic single-particle en- 
ergies is also related to the intrinsic excitation energy, 
similar to the energy variances. As a further measure 
of excitation, we introduce occupation probabilities, i. e., 
the probabilities to find a given adiabatic state within the 
space of occupied TDHF states. Adiabatic states and oc- 
cupation probabilities serve as analyzing tools, e. g., to 
investigate the nuclear Landau-Zener effect within self- 
consistent mean-field models. The scheme has been im- 
plemented numerically in fully three-dimensional TDHF 
without any symmetry restrictions and with all the terms 
of the Skyrme energy functional included. 

Two test cases of head-on collisions were considered, 
fusion of the 16 0+ 16 system at low scattering energy 
and deep inelastic scattering of 4 He+ 16 0. The newly de- 
fined boost-invariant and adiabatic single-particle ener- 
gies show the expected behaviors. For the symmetric 
16 0+ 16 system, the splitting of the asymptotically de- 
generate levels in the interaction regime is clearly seen. In 
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FIG. 3: Time evolution of the central collision of 4 He+ le O with cm. energies of 75 MeV and 125 MeV. (a) rms radius for 
the colliding system as a function of time; (b) adiabatic occupation probability P; of the last occupied neutron state; and (c) 
interacting single-particle adiabatic states. The labels of the single-particle states at t=0 are near the left axis. 



both cases, one finds the mutually complementing signals 
for intrinsic excitation and, in particular, several nicely 
developed level crossings with, in the case of 4 Hc+ 16 0, all 
signatures of a nuclear Landau- Zcncr effect. The trend 
from more adiabatic evolution at low energies to clean 
diabatic transitions at high collisional energy, c. g., is 
clearly apparent. 

These first results are very encouraging. The boost- 
invariant Hamiltonian with its single-particle energies 
and the corresponding adiabatic basis are promising tools 
for analyzing TDHF simulations of heavy-ion reactions 
and understanding their relation to the other widely used 



time-dependent method based on single-particle orbitals, 
the expansion in the adiabatic basis. 
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